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METHOD FOR RADIOLOGICAL EXAMINATION OF AN OBJECT 

DESCRIPTION 

This invention relates to a method for multiple-energy x-ray examination of an 

object. 

Radiological processes consist of having an x-ray pass through an object to be 
analyzed in order to deSHee^tbe distribution of various categories of materials in this 
object, absorbing the x-ray in different ways. A very common application is 
osteodensitometry, where the mass and density of bone tissue in a patient are 
analyzed, distinguishing these tissues from soft tissues. 

It is typical to use a broad spectrum of x-rays and to divide it into bands 
measured separately by respective measuring channels. Since the coefficients factors 
of absorption or attenuation of the x-ray by a given tissue category differ for each of 
the bands, the theoretical problem comes down to solving a system of equations, the 
number of which is the same as that of the measuring bands and each comprising two 
unknowns (the thicknesses and the masses of soft and bony tissue passed through). 
The problem becomes possible to solve after calibrating, achieved by having the x-ray 
pass through various standards provided with known thicknesses of materials with 
absorptive properties similar to those of the measuring material, in particular 
plexiglass and hydroxyapatite to simulate the soft tissue and the bony tissue. The 
mathematical parameters of a model linking the attenuation measurements to the 
thicknesses of the materials can than be calculated. 

It should be noted that in actual fact, the human body contains three main 
categories of tissue: bony tissue, lean tissue and fatty tissue, but only two of them are 
generally taken into consideration due to the difficulty of distinguishing these three 
categories in the measurements, such that lean and fatty tissue are deliberately taken 
into consideration together. Other methods are then applied to distinguish their 
proportions in the soft tissue. 

The breadth of the spectrum makes it possible to have a must larger number of 
measuring bands than would be needed to obtain a solution, and to use all of them to 
obtain more precise results by making use of all of the absorption data obtained. In the 
article "Measurement of bone material using a multiple-energy x-ray absorptiometry," 
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by J. Swanpalmer, R. Kullenberg, T. Hansson, Phys. Med. Biol., Vol. 43, 1997: pp. 
379-387, where there are 23 measuring bands and where the three categories of tissue 
are taken into consideration, it is proposed to combine the measuring groups in threes 
in every way possible to obtain 1771 (23 x 22 x 21 / 2 x 3) systems of three equations 
with three unknowns that yield as many groups of results. A criterion of choice should 
then be applied. The authors advise choosing as the true result the one that is at the 
mean value, or the median value, for the most important parameter, which can be the 
bony mass passed through. 

Contrary to the above recommendation, the authors take into consideration the 
three categories of tissue. That does not call into question the validity of their method, 
provided the patient is exposed to a much stronger radiation intensity to reduce the 
uncertainties of the measurements to acceptable proportions. One could also simply 
modify their method and apply it to measurements on only two tissues, or more 
generally speaking, two categories of materials. However, this method is flawed in 
that it does not take into account the noise on the measurements, and it thus yields 
results with noise. 

The purpose of the invention is to perfect such methods by combining 
numerous results, and it comprises improving the result selection criteria. More 
precisely, in its most general form it relates to a method for radiological examination 
of an object in which at least two categories of materials are taken into consideration, 
comprising: the use of broad spectrum x-rays; measurements of the x-rays by bands of 
the spectrum; expressions (M?)of thicknesses or masses of the two categories of 
materials passed through by the x-rays, the expressions (M?) being functions of at 
least two of the measurements (mesk) and coefficient (A); and applying a selection 
criterion from among the expressions (M?) to deduce from it an expression (final M?) 
considered true; characterized in that the selection criterion comprises a combination 
(f) of the expressions with weighting factors (a), and a calculation of weighting 
factors such that the combination has minimal noise (minimal variance in 
mathematical language) calculated according to the noise levels on the measurements 
(variance on the measurements). 

The invention will now be described referring to the figures: 

- Fig. 1 is a spectrum view, 

- Fig. 2 is a distribution of results, 
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- and Fig. 3 is a flowchart summarizing the invention, which can be referred to 
throughout the following description. 

The attenuation of the x-rays can be expressed as a function of the thicknesses 
passed through for each of the materials of indicators x and y, or of their masses M 
(density per surface unit) in the direction of the x-rays. The spectrum of 
measurements in Fig. 1 is divided into N bands generally marked by indicators i and j. 
The attenuations will vary in each of the bands according to variable absorption 
coefficients for the two materials. If we call mes; or mesj the measurements for a band 
of energy i or j, the masses passed through M x and M y can each also be expressed by 
the general formula 

(M?) = Ai + A2 • mesj + A 3 • mesj + A4 * mes 2 , -1- A5 • mes 2 , + A<$ • mesj • mesj 
The measurements considered in this example being attenuation measurements, we 
will have for each measurement channel i (corresponding to a band of the spectrum ) 
the relation mesi = In (noi/ni) where noi is the number of photons arriving on the 
object and ni the number of photons having passed through the object. Since the non- 
linearity of the functions M x and M y as a function of the measurements is slight in 
practical reality, one can make use of this second degree polynomial function that 
comprises six coefficients Ai through A 6 . The degree of the polynomial can be 
adjusted depending on the problem. For example, to analyze objects made of 
materials having atomic numbers higher than that of biological tissue, as in non- 
destructive monitoring for examining a metal object, the measurements considered in 
this example being attenuation measurements, we will have for each measurement 
channel i (corresponding to a band of the spectrum ) the relation mesi = In (noi/ni) 
where noi is the number of photons arriving on the object and ni the number of 
photons having passed through the object. 

These coefficients can be found in a phase of calibration through standards , 
sometimes called phantoms or wedges in the prior art, consisting of shaped parts with 
known, differing thicknesses made of materials simulating with their attenuation 
properties the materials of the object that will actually be measured. Each of these 
standards is thus radiated with the x-rays for a long period, making it possible to 
reduce the noise's influence on the measurements. The spectrum measured for each of 
the standards still yields N measurements resulting from the decomposition of the 
spectrum into as many bands. By then combining two series of measurements i and j 
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taken for two bands N for each of the standards, the coefficients A are analyzed to 
adjust the functions M x and M y to the measurements. In this case, where there are six 
coefficients A for each of the two functions, and where the combinations of two bands 
of measurements are being considered, the measurements for each of the bands will 
have to be for six standards in order to produce a unique solution. More numerous 
standards may also be used to improve precision on the solution. An error function 
minimizing will then be applied. 

This determining of coefficients A is repeated for various combinations of 
pairs of measurements. This was still the case in the above-mentioned previous 
article; however, it was noted that it was pointless to do all of the combinations, 
numbering Nx(N-l) I 2, to fully utilize the measurements and that (N-l) combinations 
were actually enough to gather all of the data. 

A preferred method consists of selling in the beginning the band of 
measurements with the least noise (for example the one that has the greatest signal n - 
number of photons at reception) and successively combining it with each of the other 
bands of measurements for the combinations. One finally obtains (N-l) estimates of 
the two functions M x and M y , which are noted Mi? , M 2 ? , M?N-\ for each of 
these two functions. 

At this stage of the process, the functions M x andM y representing the 
lengths passed through in two materials representing bony tissue and soft tissue can 
be converted to functions M u , M v and M w representing equivalent lengths passed 
through in the bony tissue, lean tissue and fatty tissue by combining linearly M x and 
My in three different ways determined by experiment. This conversion method is 
independent of the invention and already ready known in the prior art. 

Since there is no reason to prefer either of these M? estimates, a selection 
criterion must be applied to obtain the final M? estimate that will be considered true. 
In the previous article, one of the expressions obtained was directly selected 
according to a criterion of classification (the mean result) or median of the values 
considered by the expressions for one of the results. In the invention, the M? 
expressions will be combined; e.g. linearly according to the final M? formula = (al 
M?\) + (a2 M?2) +...+ (aN-1 Af ?N -1) while minimizing the noise; the coefficients 
ai, etc. have a sum equal to the unit 
(al +a2 + ... +aN-l = 1). 
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For each measurement channel, the noise on the number of photons follows a 
Poisson statistical law, the result of which is independent for each of the channels. 
The co variance matrix of the N-l results can be expressed according to the formula 

oM?! cM? f 1 

r,- 2 

*=i dmesk dmesk Nk 

The variance on the linear combination yielding final M? is expressed by the 
formula f = (ai, aw-i) • T • l (ai, a N -i); this quantity f attains an optimal value 
when its derivative according to all of the variables is null, that is, the noise's 
influence is minimized, i.e., 

3f / dsn = 0 



5f/5a 2 = 0 



df/da N _, = 0 

The numerical resolution of this system provides the coefficients al, a2, etc. N 
final M?, that is, the masses passed through in the two categories of materials. 

These are the operations completed in the invention; it should be noted that the 
measurements (mes x ) used in the formula yielding T y and so forth are the 
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measurements carried out through the actual object to be analyzed, but not those that 
were made by calibration to determine the coefficients A of the functions M?. 

The invention makes it possible to avoid the dilemma of prior methods in 
which part of the measurement energy was given up or, on the contrary, 
measurements that barely differed were accepted: it indeed uses the entire spectrum, 
but divides it into bands that are numerous enough to enable the measurement of each 
one to be usefully compared to other measurements made on distant bands. It thereby 
makes it possible, among other advantages, to analyze lean organisms as well as fatty 
ones. 

Lastly, the invention can be broadened to a greater number of materials than 
two, which can be advantageous in particular for imaging processes with contrast 
products, where three variables must be considered or for checking baggage 
(searching for explosives); and it can be applied when considering combinations of 
the M ? functions other than linear ones, 

Applications of the invention are: 

- osteodensitometry 

- obtaining the density of bone mass 

- obtaining body composition (fat mass, lean mass) 

- food-processing inspection, e.g. detecting bone fragments in meat or 
detecting pieces of glass in meals being served, 

- baggage inspection; searching for explosives, illegal products (weapons, 
food, drugs,...). 
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[graphs] 



RESULT AT MOYEN = MEAN RESULT 
RESULT AT MEDIAN = MEDIAN RESULT 
NUMERO DU RESULTAT - RESULT NUMBER 
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2/2 

MEASUREMENTS THROUGH The STANDARDS 
FOR A SPECTRUM WITH N BANDS (i j) 

CALCULATION OF THE COEFFICIENTS (A) 
OF THE FUNCTIONS YIELDING THE THICKNESSES 
PASSED THROUGH (KNOWN) ACCORDING TO 
THE MEASUREMENTS (mesj, me Sj ) 

MEASUREMENTS THROUGH THE OBJECT TO BE ANALYZED 



ESTIMATES (M?) OF THE THICKNESSES 
PASSED THROUGH BY APPLYING 
THE CALIBRATION STANDARDS 

COMBINATION (f) OF THE ESTIMATES 
WITH WEIGHTING FACTORS (a) 



ADJUSTMENT OF THE WEIGHTING FACTORS 
TO MINIMIZE A VARIATION OF COMBINATION 
AS A FUNCTION OF THE MEASUREMENTS 

FINAL ESTIMATE OF THE THICKNESSES 
BY APPLYING THE COMBINATION (f) 
AND THE ADJUSTED FACTORS 



FIG. 3 



